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PREFACE 


The  US  Army  Engineer  District,  Norfolk  (CENAO) ,  requested  the  Coastal 
Engineering  Research  Center  (CERC) ,  US  Army  Engineer  Waterways  Experiment 
Station  (WES) ,  to  assist  in  the  design  of  a  Beach  Erosion  and  Hurricane  Pro¬ 
tection  Project  for  Virginia  Beach,  Virginia.  The  study  was  divided  into  two 
major  parts  consisting  of  a  seawall  design  and  beach  nourishment  design.  This 
report  is  the  third  in  a  series  of  three  and  addresses  the  evaluation  and 
design  of  the  beach  fill.  Funding  authorizations  by  CENAO  were  granted  in 
accordance  with  Intra-Army  Order  No.  AD-86-3018. 

This  study  was  conducted  at  CERC  tinder  the  general  direction  of 
Dr.  James  R.  Houston,  Chief,  CERC;  Mr.  Charles  C.  Calhoun,  Jr. ,  Assistant 
Chief,  CERC;  Mr.  Thomas  W.  Richardson,  Chief,  Engineering  Development  Divi¬ 
sion;  Mr.  H.  Lee  Butler,  Chief,  Research  Division;  Mr.  Claude  E.  Chatham,  Jr., 
Chief,  Wave  Dynamics  Division;  Ms.  Joan  Pope,  Chief,  Coastal  Structures  and 
Evaluation  Branch;  and  Mr.  D.  D.  Davidson,  Chief,  Wave  Research  Branch.  This 
report  was  prepared  by  Mr.  Mark  Hansen  and  Dr.  Norman  Sheffner.  Report  edit¬ 
ing  was  performed  by  Ms.  Lee  T.  Byrne,  Information  Technology  Laboratory. 

This  study  was  closely  coordinated  with  Mr.  Dave  Pezza,  CENAO  Project 
Manager;  Mr.  Jerry  Swean,  District  Geologist;  and  Mr.  Paul  Bowen,  CENAO. 
Acknowledgment  is  made  to  all  others  involved  at  CENAO  for  their  assistance  in 
the  study. 

Commander  and  Director  of  WES  during  the  publication  of  this  report  was 
COL  Larry  B.  Fulton,  EN.  Technical  Director  was  Dr.  Robert  W.  Whalin. 
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CONVERSION  FACTORS,  NON- SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non- SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


_ Multiply _ 

cubic  yards 
feet 

miles  (US  statute) 
pounds  (mass) 


_ Bx _ 

0.7645549 

0.3048 

1.609347 

0.4535924 


To  Obtain 
cubic  metres 
metres 
kilometres 
kilograms 
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COASTAL  ENGINEERING  STUDIES  IN  SUPPORT  OF  VIRGINIA  BEACH.  VIRGINIi. 
BEACH  EROSION  CONTROL  AND  HURRICANE  PROTECTION  PROJECT 
EVALUATION  AND  DESIGN  OF  BEACH  FILL 


PART  I; 


INTRODUCTION 


Project  Description 

1.  The  proposed  Virginia  Beach,  Virginia,  Beach  Erosion  Control  and 
Hurricane  Protection  Project  is  one  of  the  largest  and  most  complex  coastal 
erosion  and  flooding  projects  in  recent  US  Army  Corps  of  Engineers  experience. 
The  City  of  Virginia  Beach  is  located  on  the  east  coast  of  the  United  States 
just  south  of  the  entrance  to  Chesapeake  Bay  (Figure  1).  The  project  area 
consists  of  6  miles*  of  heavily  developed  commercial  and  urban  shoreline  that 
extends  north  from  Rudee  Inlet  to  89th  Street  (Figure  2) .  This  shoreline  is 
subject  to  severe  damages  from  both  hurricanes  and  extreme  extratropical 
storms.  Both  the  August  1933  hurricane  and  the  March  1962  extratropical  storm 
("the  Ash  Wednesday  stoim")  devastated  this  coastal  area.  Storm  damages  have 
included  loss  of  the  beach,  destruction  of  the  bulkhead  and  seawall  system, 
damage  to  buildings,  and  inshore  flooding.  In  addition,  beach  erosion  has 
been  a  continuing  problem.  Since  1962,  annual  harbor  dredging  and  pumping 
operations  to  by-pass  sand  at  Rudee  Inlet  and/or  the  trucking  in  of  sand  from 
other  sources  have  been  sponsored  by  the  Federal,  State,  and  city  governments 
to  maintain  a  beach  width  of  approximately  65  ft  and  a  crest  elevation  of 
+5.4  ft.** 

2.  Existing  protection  consists  of  a  combination  of  various  bulkheads 
with  crest  elevations  between  10  and  12  ft  NGVD  and  nourished  beach.  In  1970 
the  US  Army  Engineer  District,  Norfolk  (CENAO) ,  completed  a  feasibility  study 
that  recommended  construction  of  a  sheet-pile  seawall  with  a  concrete  cap  at 
elevation  15  ft  NGVD  and  heavy  stone  at  the  base.  By  1983,  results  of  the 
previous  study  had  been  reevaluated  and  incorporated  into  an  initial  (Phase  I) 
seawall  design  and  beach  erosion  concept.  The  seawall  was  designed  with 


*  A  table  of  factors  for  converting  non- SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 

**  All  elevations  cited  herein  are  in  feet  referred  to  the  National  Geodetic 
Vertical  Datum  (NGVD) . 
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Figure  1.  Location  map  of  study  area 


guid*>-  -e  from  the  Shore  Protection  Manual  (SPM)  (1984)  based  primarily  on 
monochromatic  wave  theory.  Adequate  storm  protection  was  to  be  provided  by 
the  seawall  without  sacrificing  aesthetics  of  the  ocean  view. 

3.  The  proposed  construction  project  is  a  new  stepped-face  seawall  with 
curved  parapet,  located  just  seaward  of  the  existing  seawall  (between  Rudee 
Inlet  and  57th  Street) .  The  existing  dune  field  will  be  raised  and  widened 
from  57th  Street,  north  to  89th  Street.  Both  flood-control  structures  will  be 
fronted  by  a  continuously  maintained  beach  berm. 

Study  Background 

4.  This  report  is  the  last  of  a  series  of  three  reports  on  Coastal 
Engineering  Studies  which  were  conducted  by  the  Coastal  Engineering  Research 
Center  (CERC) ,  US  Army  Engineer  Waterways  Experiment  Station,  to  assist  the 
CENAO  in  the  Advanced  Engineering  and  Design  of  the  Virginia  Beach,  Virginia, 
Beach  Erosion  Control  and  Hurricane  Protection  Project.  The  other  two  reports 
(Heimbaugh  et  al .  1988;  Lillycrop,  Pope,  and  Abel  1988)  consist  of  the  seawall 
design  (i.e.,  physical  model,  overtopping  tests,  and  physical  model  pressure 
or  wave  loading  tests)  and  selection  of  the  design  parameters  and  results  of 
the  physical  model  seawall  overtopping  tests.  This  report  evaluates  the 
effectiveness  of  various  beach- fill  designs  and  concludes  with  design  of  the 
selected  plan.  Figure  3  presents  a  flowchart  of  the  coastal  engineering 
studies . 

5.  Selection  of  design  waves,  storm  surge  hydrographs,  and  runup¬ 
overtopping  rates  was  crucial  to  the  development  of  the  most  hydraulically 
efficient  seawall  geometry  and  in  defining  the  short-term  beach  stability. 
Coastal  engineering  studies  consisted  of  selection  of  design  storms  from  the 
historical  record,  simulating  the  wave  field  for  each  of  these  storms,  estab¬ 
lishing  the  design  surge  hydrographs,  and  developing  a  two-dimensional  (2-D) 
hydrographic  model  to  test  and  measure  overtopping  rates. 

6.  The  seawall  overtopping  tests  were  concerned  with  hydraulically 
designing  the  most  efficient  seawall  plan  and  developing  overtopping  rates 
that  could  be  used  in  interior  flooding  design.  The  study  used  a  2-D  physical 
model  (Heimbaugh  et  al.  1988)  over  a  range  of  sea  state  conditions  to  measure 
wave-induced  overtopping.  Parameters  which  required  evaluation  and  were 
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Figure  3 .  Flowchart  for  coastal  engineering  studies 
Virginia  Beach,  Virginia 


incorporated  into  the  model  test  were  design  wave  conditions,  the  design  storm 
surge  levels,  seawall  geometry,  and  beach  profile. 

7.  This  report  is  divided  into  two  parts.  The  first  part  discusses 
results  of  numerical  modeling  using  a  profile  response  model  to  evaluate  the 
effectiveness  of  various  beach* fill  designs.  Based  upon  results  of  the  numer¬ 
ical  modeling  effort,  standard  empirical  procedures,  and  field  data,  the  sec¬ 


ond  part  contains  engineering  parameters  of  the  recommended  beach-fill  design. 
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PART  II.  BEACH-FILL  DESIGN  EVALUATION 


8.  The  objective  of  this  section  is  to  determine  if  the  beach- fill  and 

dune  designs  proposed  by  the  CENAO  are  sufficient  for  protecting  the  low-lying 
backshore  areas  from  storm -induced  flooding.  The  area  of  concern  extends  from 
Rudee  Inlet  to  88th  Street  of  Virginia  Beach,  Virginia.  Two  designs  are  eval¬ 
uated:  one  with  a  protective  seawall  (Figure  4)  and  one  without  a  seawall 

(Figure  5).  In  each  case,  a  100-ft  flat  berm  is  specified  in  the  design 
profile. 

9.  The  numerical  model  used  to  evaluate  the  proposed  design  is  the 
modified  Kriebel  beach  erosion  model.  The  basic  theory  of  the  model  is 
described  in  detail  by  Kriebel  (1984a,  b) ,  who  provides  a  program  listing. 
Birkemeier  et  al.  (1987)  discuss  model  strengths  and  limitations.  A  modifica¬ 
tion  was  made  at  CERC  to  expand  the  capability  of  the  model  to  Include  the 
simulation  of  vertical  seawalls.  This  enhancement  will  be  subsequently  dis¬ 
cussed.  A  first  fundamental  assumption  of  the  model  is  that  the  offshore 
depth  can  be  described  by  an  equilibrium  relationship  of  the  following  form: 

h(x)-Ax^/^  (1) 

where 

h  -  depth  below  mean  sea  level  (MSL) 

A  -  equilibrium  coefficient 

X  -  distance  offshore 

This  relationship  was  first  reported  by  Bruun  (1954)  and  was  substantiated  by 
Dean  (1977)  in  a  study  in  which  Equation  1  was  shown  to  satisfactorily 
describe  502  offshore  profiles  measured  along  the  eastern  coast  of  the  United 
States.  Subsequent  evaluations  by  Hughes  (1978)  and  Moore  (1982)  have  shown 
the  equation  to  be  valid  for  beach  profiles  in  the  laboratory  and  for  beaches 
in  various  parts  of  the  world. 

10.  The  above  studies  have  shown  that  the  equilibrivim  coefficient  A 
is  correlated  to  the  mean  grain  size  of  the  beach  material.  For  example,  a 
grain  size  of  0.25  mm  corresponds  to  an  A  value  of  approximately  0.13  m^^® 
(0.19  ft^^®) .  Figure  6,  reproduced  from  Moore  (1982),  shows  the  relationship 
of  grain  diameter  to  equilibrium  coefficient.  The  preferred  procedure  for 
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57  th  TO  69  th  STREETS 

Figure  5.  Typical  design  cross  section  without  protective  seawall 

with  100-ft-wide  berm 
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Sediment  Diameter  (nm) 

Figure  6.  A  versus  sediment  diameter  (after  Moore  1982) 

determining  an  A  value  for  a  given  location  Is  to  use  the  existing  offshore 
bathymetry  to  compute  an  average,  or  best  fit,  coefficient.  If  the  existing 
profile  cannot  be  considered  to  be  In  equilibrium  (l.e.,  a  storm  event  has 
recently  occurred  or  coastal  structures  are  In  the  Immediate  vicinity)  or  If 
no  profile  data  are  available,  the  mean  grain  size  (D^g)  can  be  used  to  esti¬ 
mate  a  value  for  the  coefficient  A  according  to  Figure  6. 

11.  The  equilibrium  profile  of  the  form  of  Equation  1  was  shown  by  Dean 
(1977)  to  result  from  the  assumption  that  the  equilibrium  profile  Is  a  func¬ 
tion  of  the  dissipation  of  wave  energy  per  unit  volume  of  water  in  the  active 
surf  zone.  Transport  of  sediment  In  the  surf  zone  Is  assumed  to  be  a  primary 
function  of  the  dissipation  of  the  flux  of  this  wave  energy.  This  basic 
assumption  Is  used  to  formulate  a  sediment  transport  equation  In  the  following 
form: 


Q,  -  k(D  -  D„)  (2) 

where 

0  -  energy  dissipation  written  as 
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D 


1  a (EC,) 
h  dx 


E  -  total  wave  energy 
C,  -  group  velocity 

The  parameter  D,,  represents  the  dissipation  for  an  equilibrium  profile  in 
which  the  depth  and  offshore  distance  are  related  according  to  Equation  1. 

The  parameter  k  is  an  empirical  coefficient  found  by  Moore  (1982)  to  have  a 
value  of  2.2  X  10"®  m'/N  (0.001144  ftVlb)  .  Linear  wave  theory  defines  the 
two  dissipation  relationships  as  follows; 

°  “  it 

°«q  “24 


In  Equations  3  and  4,  p  represents  the  fluid  density,  g  represents  the 
acceleration  of  gravity,  and  7  represents  the  ratio  of  wave  height  H  to 
depth  h  . 

12.  The  above  derivations  assume  the  surf  zone  is  dominated  by  spilling 

breakers  such  that  the  breaking  wave  height  is  a  constant  fraction  of  the 
depth,  l.e.  H  -  yh  .  Equations  2,  3,  and  4  are  used  to  compute  a  net  volume 
of  erosion  or  deposition  for  the  area  bounded  approximately  by  the  shoreline 
(MSL)  and  the  breaker  line.  Note  that  for  a  beach  in  equilibrium,  no  trans¬ 
port  of  material  occurs,  since  D  -  .  The  basic  premise  of  the  model  is 

that  this  volume  of  material  is  either  supplied  from  the  dune  and  berm  or 
deposited  on  the  berm  depending  on  whether  erosion  or  deposition  is  indicated 
offshore.  Deposition  on  the  dune  face  during  poststorm  recovery  is  not  per¬ 
mitted.  The  assumption  is  therefore  made  that  erosion  of  the  dune  is 
irreversible. 

13.  The  equilibrium  profile  concept  has  one  disadvantage;  it  indicates 
that  the  offshore  depth  monotonically  increases  with  offshore  distance.  This 
assumption  precludes  the  formation  of  offshore  bars  or  troughs  which  are  known 
to  occur  in  many  locations;  however,  the  average  profile  has  been  shown  to  be 
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adequately  described  by  this  equilibrium  profile  concept  (Dean  1977) .  The 
primary  goal  of  the  dune  erosion  model  is  not  to  compute  offshore  erosion 
patterns,  but  to  compute  the  above  MSL  changes  that  result  from  a  specific 
storm  event  based  on  the  total  volume  of  erosion  (or  deposition)  occurring 
between  the  shoreline  and  the  breaker  zone.  The  use  of  an  average  offshore 
profile  to  compute  this  total  volume  of  erosion  is  therefore  a  reasonable 
assumption  for  areas  which  are  not  dominated  by  offshore  anomalies  that  pro¬ 
duce  highly  variable  longshore  transport  patterns. 

14.  The  present  equilibrium  profile  approach  to  dune  erosion  modeling 
also  asstmes  that  alongshore  transport  is  constant  and  that  erosion  of  the 
dune  face  is  solely  a  function  of  the  onshore -offshore  transport.  This  one¬ 
dimensional  (1-D)  (onshore-offshore)  approximation  of  a  2-D  (onshore- 
offshore/alongshore)  natural  process  is  an  acceptable  one  when  alongshore 
transport  is  in  balance.  For  example,  sediment  entering  the  test  section  is 
assumed  to  leave  the  test  section  in  the  identical  spatial  and  temporal  dis¬ 
tribution  such  that  the  net  alongshore  effect  is  negligible  on  the  profile. 
This  assumption  is  a  reasonable  one  for  the  Virginia  Beach  project,  which  is 
characterized  by  a  relatively  straight  coastline  free  of  major  structures  such 
as  jetties  or  breakwaters.  Although  the  1-D  simplification  is  well  suited  for 
the  Virginia  Beach  area,  the  potential  impact  of  the  approximation  on  the 
computed  results  must  be  critically  evaluated.  For  example,  breaching  of  the 
dune  complex  was  experienced  during  the  1962  northeaster.  This  undoubtedly 
occurred  because  of  variations  in  the  dune  height  along  the  dune  crest,  such 
as  walk-over  areas  and  locations  of  channelized  wind  erosion,  in  which  the  2-D 
assumption  is  violated.  Similarly,  the  net  alongshore  littoral  drift  has  been 
reported  to  be  to  the  north.  The  extreme  southern  limit  of  beach  fill  could 
become  sediment  starved  as  the  newly  placed  material  moves  north  with  no  addi¬ 
tional  material  to  replace  it  from  the  south.  In  this  case,  that  1-D  assump¬ 
tion  would  also  be  violated  since  the  alongshore  inflow  would  not  equal  the 
alongshore  outflow.  These  types  of  analyses  must  be  considered  in  the  evalua¬ 
tion  of  the  results  of  any  idealized  numerical  model.  As  with  all  numerical 
models,  it  is  the  responsibility  of  the  user  to  correctly  interpret  the  numer¬ 
ical  results.  This  is  especially  true  for  sediment  transport  models  in  which 
the  governing  physics  are  not  well  understood,  and  empirical  relationships  are 
used  to  describe  the  transport  process.  This  approach  is  acceptable  as  long 
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as  the  basic  asstmptlons  used  in  the  determination  of  the  relationships  are 
not  severely  violated. 

15.  A  quantitative  evaluation  of  the  dune  erosion  model  was  made  by 
Birkemeier  et  al.  (1987),  in  which  14  prestorm  and  poststorm  surveys  repre¬ 
senting  four  separate  storm  events  were  used  to  test  the  capability  of  the 
model  to  realistically  predict  volumes  of  erosion  resulting  from  known  storm 
surges.  Results  of  the  comparison  showed  that  the  average  percentage  of  devi¬ 
ation  between  calculated  erosion  volumes  above  MSL  and  measured  volumes  was 
109  percent.  This  value  represents  the  average  of  five  underpredictions  aver¬ 
aging  55  percent  of  the  measured  and  eight  oveirpredictlons  representing 
145  percent  of  measured.  One  profile  of  the  14  was  influenced  by  adjacent 
structures  and  was  not  used  in  the  above  percentages.  The  profiles  used  for 
this  analysis  did  not  Include  profiles  with  seawalls;  therefore,  an  additional 
evaluation  of  the  model  was  made  by  Kraus  et  al.  (in  preparation).  Measured 
elevation  changes  in  front  of  a  seawall  resulting  from  a  documented  storm 
surge  were  compared  with  elevation  changes  as  computed  by  the  model  with  the 
seawall  modification.  Three  prestorm  and  poststorm  surveys  were  used  to  test 
the  model.  Average/  predicted  elevation  changes  for  the  three  profiles  were 
-0.31/-0.26,  -0.50/-0.42,  and  -0.32/-0.04  m,  respectively.  Total  volumes  of 
erosion  could  not  be  compared  since  adequate  prestorm  and  poststorm  data  were 
not  available.  Conclusions  of  both  verification  comparisons  of  the  model  show 
it  to  be  capable  of  acceptably  predicting  storm-related  erosion.  In  view  of 
the  fact  that  a  natural  variability  in  observed  erosion  occurs  on  what  would 
appear  to  be  a  straight  and  homogeneous  section  of  beach,  the  predictions  of 
the  model  to  be  approximately  ±50  percent  of  the  observed  values  represent  a 
positive  feature  of  the  model.  A  worst  case  scenario  can  be  made  by  applying 
a  "variability  factor"  of  2.0  to  the  model  predictions.  Support  of  a  factor 
of  this  magnitude  was  reported  by  Birkemeier  et  al.  (1987)  and  Chiu  and  Dean 
(1986).  The  dune  erosion  model  has  been  shown  to  effectively  predict 
quantities  of  erosion  resulting  from  single  storm  events  of  known  surge  level 
and  duration.  The  initial  analyses  that  follow  are  also  made  for  single  storm 
events  in  which  the  initial  beach  profile  is  assumed  to  be  the  design  profile. 
If  multiple  storm  events  occur  between  beach  renourishment,  a  corresponding 
increase  in  erosion  of  the  dune  should  be  anticipated.  A  limited  analysis  of 
multiple  events  is  made  to  address  this  additional  volume. 
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16.  Detailed  profiles  from  19  locations  along  Virginia  Beach  were  used 
to  compute  an  average  equlllbrlvim  coefficient  for  use  In  the  model.  The  Indi¬ 
vidual  values  of  the  coefficients  for  each  of  the  profiles  are  shown  In 
Table  1.  As  can  be  seen,  the  coefficients  Indicate  that  all  profiles  are  very 
similar  In  shape.  This  result  demonstrates  that  the  project  area  Is  well 
suited  to  the  1-D  analysis  of  the  dune  erosion  model.  All  of  the  above  pro¬ 
files  were  supplied  to  CERC  by  CENAO  and  represent  profiles  collected  on 
20  June  1986.  The  computed  range  of  coefficient  values  Is  0.165  to 
0.190  ft^^^,  corresponding  to  a  sediment  diameter  of  0.25  to  0.30  mm.  An 
approximate  average  value  of  0.18  ft^^^  was  used  for  all  computations.  This 
value  corresponds  to  a  mean  sediment  diameter  of  0.28  mm,  which  Is  consistent 
with  the  sediment  analysis  provided  CENAO  by  Waterways  Surveys  &  Engineering 
Ltd.  (1984)  for  the  Cape  Henry  area.  Material  of  this  diameter  Is  also  con¬ 
sistent  with  the  following  potential  sources  of  fill  material  reported  by 
CENAO  (1984):  Llnkhorn  Bay  (0.29  mm).  The  Narrows  and  Broad  Bay  (0.25  to 
0.34  mm).  Fort  Story  and  Seashore  State  Park  (0.26  nun),  Lynnhaven  Inlet 
(0.28  mm).  Ocean  Naval  Air  Station  area  (0.34  mm  north,  0.26  mm  south).  Char¬ 
ity  and  Pungo  Ridges  (0.32  mm).  Coarser  material  was  reported  for  Bonney's 
Comer  (0.39  mm)  and  the  east  end  of  the  Thimble  Shoal  channel  (0.42  mm), 
while  finer  material  was  Indicated  offshore  of  the  Virginia  Beach  project 
(0.19  mm)  and  the  Dam  Neck  disposal  site  (0.16  mm). 

Table  1 

Equilibrium  Coefficients 


Coeff.  Coeff.  Coeff. 


Street 

Street 

ftV3 

Street 

ftV3 

2 

0.165 

5 

0.165 

8 

0.170 

11 

0.180 

14 

0.180 

16 

0.180 

18 

0.170 

21 

0.175 

24 

0.180 

27 

0.180 

30 

0.190 

33 

0.180 

35 

0.180 

38 

0.180 

40 

0.185 

43 

0.180 

46 

0.180 

49 

0.185 

52 

0.180 
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17.  The  two  beach- fill  designs  provided  by  CENAO  were  schematized  to 
best  represent  the  design  profiles  shown  in  Figures  4  and  5.  These  data  were 
input  to  the  dune  erosion  model  for  all  simulations.  The  NGVD  and  HSL  are 
equal  in  value  for  the  project  area,  and  all  elevations  are  referenced  to 
these  datum.* 

18.  Four  storm  events  of  record  were  supplied  by  CENAO  for  evaluation 
of  the  two  proposed  beach- fill  designs.  A  request  was  also  made  by  CENAO  to 
test  the  design  dune  system  with  a  hurricane  having  a  peak  surge  elevation  of 
9.5  ft.  This  test  storm  was  assembled  by  scaling  the  1933  hurricane  such  that 
the  peak  surge  was  linearly  scaled  up  from  8.7  to  9.5  ft.  The  approximate 
return  period  for  each  storm  event  was  extracted  from  the  adopted  Virginia 
Beach  elevation-exceedance  curves  published  by  CENAO  (1983)  and  shown  in  Fig¬ 
ure  7.  Pertinent  storm-related  data,  representing  three  northeasters  and  two 
hurricanes,  are  shown  in  Table  2.  The  surge  profiles  are  shown  in  Appendix  A. 

19.  Each  numerical  simulation  generates  a  complete  updated  cross  sec¬ 
tion  at  each  computational  time  step.  To  summarize  these  results  in  a  concise 
and  meaningful  manner,  two  output  values  are  reported  for  each  simulation. 
First,  the  maximum  computed  poststorm  recession  of  the  dune/berm  face  is  pro¬ 
vided  as  a  meaningful  indicator  of  potential  storm  damage  by  structural  under¬ 
mining.  The  second  output  represents  the  net  volume  of  erosion  per  foot  of 
beach  width  (alongshore)  that  results  from  the  total  stoirm  event.  This  figure 
reflects  the  fact  that  beach  recovery  does  occur  and  indicates  the  amount  of 
beach  renourishment  which  would  be  required  to  restore  the  beach  to  prestorm 
conditions.  Results  of  the  initial  simulations  are  shown  in  Table  3. 

20.  The  prestorm  and  poststorm  profiles  for  each  of  the  above  storm 
events  for  both  the  with  and  without  seawall  designs  are  shown  in  Appendix  B. 
Note  that  for  the  storm  events  with  surges  greater  than  5.5  ft,  no  erosion  of 
the  dune  face  is  indicated  in  spite  of  the  submergence  of  the  berm.  This  is  a 
consequence  of  the  numerical  algorithm  in  which  the  assumption  is  made  that 
the  flat  berm  section  erodes  first  before  the  dune  erodes.  If  the  surge  level 
remains  in  contact  with  the  dune  face  for  a  substantial  portion  of  the  storm 
event  without  completely  eroding  the  flat  berm,  this  assumption  may  represent 
a  poor  approximation  of  the  natural  process.  In  this  case,  an  evaluation  of 
the  results  should  be  made  as  mentioned  above.  Note  also  in  Table  3  that  a 


*  Personal  Communication,  1987,  R.  Owen  Reece,  Jr.,  CENAO,  Norfolk,  VA. 
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EXCEEDENCE  FREQUENCY  IN  PERCENT 

Figure  7.  Exceedance  frequency  curves  for  Virginia  Beach,  Virginia 


Table  2 


Design  Storms 


Peak  Surge 

Duration 

Return  Period* 

Tvne 

ft 

hr 

vears 

August  1933 

hurricane 

8.7 

23.0 

100 

Test  storm 

hurricane 

9.5 

23.0 

200 

April  1956 

northeaster 

6.3 

71.0 

10 

March  1962 

northeaster 

6.7 

78.0 

15 

April  1978 

northeaster 

5.8 

56.0 

5 

*  Based  on 

still -water  level  of 

the  surge. 

Table  3 

Volume  of  Erosion 

Maximum  Recession 

Design 

ft 

vd  /ft 

August  1933 

no  seawall 

-50.63 

7.87 

seawall 

-50.72 

7.87 

Test 

no  seawall 

-56.76 

8.50 

seawall 

-56.89 

8.50 

April  1956 

no  seawall 

-25.96 

4.97 

seawall 

-25.94 

4.98 

March  1962 

no  seawall 

-42.27 

9.45 

seawall 

-42.27 

9.45 

April  1978 

no  seawall 

-41.71 

8.27 

seawall 

-41.70 

8.27 

small  discrepancy  is  shown  in  the  with  and  without  seawall  results.  This 
difference  is  a  result  of  the  fact  that  a  flat  berm  cannot  be  specified  as 
input  to  the  model  for  the  seawall  case.  Computationally,  a  1:200  slope  is 
used  for  this  area,  which  introduces  small  discrepancies  in  value. 

21.  The  nximerical  results  shown  in  Table  3  indicate  that  the  design 
storms  do  not  cause  erosion  of  the  entire  flat  berm  region.  In  fact,  the 
9. 5- ft  test  storm  resulted  in  a  maximum  recession  of  only  approximately  57  ft. 
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As  seen  in  figures  from  Appendix  B,  a  recession  of  this  magnitude  does  not 
reach  the  base  of  the  dune  or  the  face  of  the  seawall.  As  previously  stated, 
however,  a  worst  case  scenario  can  be  surmised  by  applying  a  variability 
factor  of  2.0  to  these  maximum  recession  values.  Even  with  this  factor 
Included,  the  design  profiles  provided  by  CENAO  are  still  shown  to  provide 
adequate  protection  to  the  dune  and  seawall  since  erosion  of  the  base  on  the 
dune  or  face  of  the  seawall  will  Just  have  begun.  To  test  this  scenario,  a 
multiple  storm  simulation  was  made.  Four  test  storms  were  constructed  by 
assuming  that  a  single  storm  event  was  Immediately  followed  by  an  additional 
storm  event.  The  storms  used  were  the  9.5-ft  test  hurricane  and  the  March 
1962  northeaster.  The  four  permutations  of  these  storms  (33-33,  33-62,  62-33, 
and  62-62)  were  subjected  to  each  of  the  CENAO  designs.  Results  of  these 
simulations  are  shown  In  Table  4.  Appendix  B  represents  the  prestorm  and 
pdststorm  profiles.  In  only  two  cases  (33-33  and  62-33)  are  erosion  of  the 
crest  of  the  dune  Indicated.  These  recession  values  are  1.94  and  1.78  ft 
respectively.  In  the  respective  seawall  cases,  a  volumetric  loss  of  the  berm 
elevation  up  to  2  ft  occurred  over  Its  entire  width  to  the  base  of  the 
seawall .  Table  4 

Model  Results.  Multiple  Simulation 


Storms 

Deslan 

Maximum  Recession 
ft 

Volume  of  Erosion 
vd^/ft 

33-33 

no  seawall 

-74.68 

13.23 

seawall 

-82.75 

13.23 

33-62 

no  seawall 

-53.81 

12.67 

seawall 

-53.82 

12.67 

62-33 

no  seawall 

-75.77 

14.00 

seawall 

-83.28 

14.00 

62-62 

no  seawall 

-56.92 

13.60 

seawall 

-56.94 

13.60 

22.  An  additional  set  of  simulations  was  made  In  which  the  berm  was 
reduced  to  a  SO- ft  width  In  order  to  assess  the  computed  erosion  resulting 
from  the  9.5-ft  surge  design  storm.  Without  the  protection  of  the  seawall, 
approximately  5.0  ft  of  the  dune  crest  would  be  eroded,  with  a  maximum 
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recession  on  the  dune/berm  face  of  34.75  ft.  Again,  results  for  the  seawalled 
case  indicate  a  volumetric  loss  from  the  berm  resulting  in  an  elevation  loss 
up  to  2  ft  over  its  entire  width  to  the  base  of  the  seawall.  Two  additional 
simulations  were  made  to  evaluate  a  case  in  which  the  berm  sloped  from  the 
dune  toe  (elevation  of  5.5  ft)  to  MSL.  Similar  results  are  obtained  for  these 
cases,  erosion  of  5.0  ft  of  the  dune  face  and  crest  without  a  seawall  and 
minor  erosion  of  the  berm  to  the  face  of  the  seawall. 

23.  Additional  simulations  were  made  to  demonstrate  the  effects  of  the 
9.5-ft  test  storm  on  variations  of  the  foreshore  slope  of  the  duned  design 
profile  and  of  the  equilibrium  profile  (which  is  an  indirect  measure  of  the 
effect  of  the  grain  size).  Table  5  presents  these  data.  As  seen  in  Table  5, 
the  effects  of  changing  the  subject  parameters  are  qualitatively  predictable. 
For  example,  a  steeper  foreshore  absorbs  more  storm  energy  and  is  therefore 
subject  to  greater  erosion.  Similarly,  a  smaller  grain  size  erodes  more 
easily  than  a  coarse  grain,  as  is  indicated  in  the  table.  None  of  the  effects 
are  seen  to  be  dramatic,  a  result  which  Indicates  that  the  CENAO  design  is  a 
stable  one  with  respect  to  the  geometric  parameters  and  sediment  size. 

Table  5 

Geometrical  Effects  on  Volume  of  Erosion 


Case  of  Erosion 
(with  9.5-ft  Surge) 

Design 

Maximum  Recession 
ft 

Volume 

vd^/ft 

Slope  -  1:20 

no  seawall 

-56.76 

8.50 

Slope  -  1:25 

no  seawall 

-40.76 

7.47 

Slope  -  1:30 
with  A  0.18 

no  seawall 

-29.69 

6.67 

A  -  0.15  (D50  -  0.10  mm) 

no  seawall 

-74.26 

8.75 

A  -  0.18  (D50  -  0.28  mm) 

no  seawall 

-56.76 

8.50 

A  -  0.21  (D50  -  0.32  mm) 
with  slope  -  1:20 

no  seawall 

-43.49 

8.31 

24.  The  numerical  results  presented  in  this  report  indicate  that  the 
dune/seawall  profile  design  developed  by  CENAO  is  adequate  for  protection  of 
the  areas  shoreward  of  the  dune  line.  In  no  case  was  the  dune  crest  signifi¬ 
cantly  eroded,  nor  was  the  seawall  uncovered  even  when  multiple  design  storms 
were  assumed  to  follow  in  succession.  Use  of  the  variability  factor  still 
shows  the  CENAO  design  to  be  capable  of  withstanding  the  design  storms  with 
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little  or  no  damage  to  the  dunes  or  seawall.  Results  indicate  that  the  energy 
contained  in  the  specified  storms  is  not  adequate  for  moving  the  volume  of 
fill  material  contained  in  the  design  profiles. 


21 


PART  III:  BEACH-PILL  DESIGN 


25.  Part  II  results  Indicated  that  a  berm  width  of  100  ft  at  elevation 
+5.4  NGVD  would  provide  adequate  protection  against  extremal  storm  events. 

The  purpose  of  this  study  was  to  design  a  beach  nourishment  project  using  the 
recommended  berm  dimensions.  To  design  the  project,  this  study  statistically 
analyzed  native  and  borrow  sediment  data,  evaluated  beach  profiles  for  design 
purposes,  calculated  overfill  and  renourishment  factors,  and  computed  design 
and  construction  volumes  and  templates. 

Methodology 

26.  The  volume  of  material  required  to  build  the  beach  to  the  recom¬ 
mended  dimensions,  i.e.  100-ft  berm  at  +5.4  ft  NGVD,  was  computed  based  upon 
the  profile  translation  method  (Delft  Hydraulics  Laboratory  1987;  Hansen  and 
Lillycrop  1988).  This  method  assumes  that  hydrodynamically  stable  fill  mate¬ 
rial  placed  on  the  subaerial  beach  will  eventually  be  transported  seaward 
until  a  new  equilibrium  profile  is  established.  The  limit  of  significant 
onshore -offshore  transport  is  assumed  to  be  located  at  the  average  profile 
closure  depth.  Representative  profiles  of  the  study  area  are  required  in  this 
procedure . 

27.  Three  assumptions  are  made  for  this  method  to  be  viable.  The  first 
assumption  is  that  all  hydrodynamically  unstable  fill  material  will  be  win¬ 
nowed  and  transported  either  seaward  of  the  profile  closure  depth  or  downdrift 
from  the  project  site.  This  winnowing  process  is  assumed  to  produce  an  over¬ 
all  resultant  grain  size  equivalent  to  the  native  material.  The  second 
assumption  is  that  grain  size  generally  controls  profile  steepness  (Bascom 
1951,  Shepard  1963).  If  these  two  assumptions  hold  true,  it  can  be  assumed, 
in  most  cases,  that  the  new  equilibrium  profile  will  be  congruent  to  the 
existing  profile,  however  translated  seaward.  The  last  assumption  is  that  the 
translated  profile  will  intersect,  or  "tie-in,"  to  the  existing  profile  at  an 
average  seaward  depth  of  extreme  onshore -offshore  transport,  often  termed 
profile  closure  (Hallermeler  1981). 

28.  As  a  result  of  the  continuous  nourishment  program  at  Virginia 
Beach,  it  is  difficult  to  ascertain  representative  profiles  that  are  in 
dynamic  equilibrixjm  with  the  local  wave  climate.  Therefore,  the  approach 
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taken  for  this  project  was  to  use  the  profile  data  set  which  possessed  the 
narrowest  berm  width.  Translating  profiles  exhibiting  a  narrow  berm  width 
would  provide  the  most  conservative  (greatest)  volumetric  requirements  for  the 
project.  Four  profile  data  sets  provided  by  CENAO  to  document  the  typical 
seasonal  variability  and  short-term  trends  (March  1986,  September  1986,  March 
1987,  September  1987).  The  March  1987  data  set  was  generally  found  '■o  have 
the  narrowest  berm  width  (Appendix  C) .  Therefore,  using  the  March  1987  pro¬ 
files  in  the  profile  translation  method  may  provide  the  maximvun  benefits  in 
design  of  a  storm  protection  structure  during  extreme  storm  events. 


29.  Unique  compartments,  or  subreaches,  are  typically  identified  in 
beach  nourishment  projects  to  divide  the  beach  into  zones  that  may  exhibit 
anomalous  erosion  rates  or  profile  shape.  Subreaches  are  often  treated  sepa¬ 
rately  when  computing  design  parameters,  i.e.  design  volvune,  composite  grain 


size.  At  the  request  of  CERC,  CENAO  identified  seven  unique  subreaches  along 
the  project  length.  Each  compartment  was  represented  by  an  established  pro¬ 
file  line  which  generally  bisected  the  subreach.  Table  6  Identifies  the  seven 
profiles  and  length  of  each  subreach,  designated  by  street  number,  used  for 
designing  the  hurricane  protection  project.  The  beach  north  of  65th  Street 
has  an  average  berm  width  that  exceeds  the  design  berm  width  (100  ft)  computed 
in  Part  II  of  this  report.  Therefore,  all  design  calculations  were  computed 

Table  6 


Leneth  of  Each  Subreach  Designated  bv  Street  Location 


Profile  Street 

Length,  ft 

5  th 

2,500  ; 

11th 

5,280 

29th 

5,300 

38th 

3,625 

49  th 

1,975 

55th 

1,850 

61  St 

2,560 

Total  23,090 
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for  the  beach  area  which  extends  from  Rudee  Inlet  northward  to  1,660  ft  north 
of  61st  Street,  l.e.  65th  Street.  Total  length  of  the  beach-fill  design  was 
23,090  linear  ft. 

30.  A  determination  of  profile  closure  depth  is  required  when  applying 
the  profile  translation  method  for  beach-fill  design.  As  stated  previously, 
it  is  assumed  the  new  equilibrium  beach  profile  will  "tie-in"  to  the  existing 
profile  at  a  depth  of  profile  closure.  Hallermeler 's  (1981)  method  was  used 
to  compute  closure  depth  (d)  at  Virginia  Beach.  The  nearshore  limit  of  sedi¬ 
ment  transport  as  given  by  Hallermeler  was  determined  by: 

d  -  2.28  H  -  68.5  (5) 


where 

d  -  limit  of  intense  onshore -offshore  transport 
(Mean  Lower  Low  Water  (MLLU)) 

H,  -  extreme  significant  wave  height  exceeding  12  hr /year  defined  by: 

H,  -  H,  +  5.6a  (6) 

where 

H,  -  significant  wave  height 
a  -  standard  deviation  of  significant  wave  height 
g  -  acceleration  of  gravity 
T,  -  corresponding  wave  period 

31.  Twenty  years  of  wave  hindcast  data  from  the  Wave  Information  Study 
(WIS)  (Jensen  1983)  Phase  II,  Station  77,  were  used  in  the  computation  of  clo¬ 
sure  depth.  For  Virginia  Beach,  the  average  significant  wave  height  was 
1.8  ft,  with  a  standard  deviation  equal  to  1.87  ft  and  an  average  wave  period 
of  5.7  sec.  Substituting  Equation  6  in  Equation  5  and  applying  the  20-year 
WIS  wave  summary  data  yield  an  average  annual  profile  closure  depth  of  approx¬ 
imately  -19  ft  NGVD  (-18  ft  MLLW). 
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32.  Profile  data  for  the  March  1986,  September  1986,  March  1987,  and 
September  1987  surveys  were  digitized  and  recorded  in  the  Interactive  Survey 
Reduction  Program  (ISRP)  (Birkemeier  1984)  format.  A  computer  program  was 
developed  to  translate  the  representative  profiles  a  specified  distance  from 
the  baseline,  at  a  given  berm  height,  and  intersecting  the  existing  profile  at 
a  given  closure  depth.  As  specified  by  results  from  Part  II  of  this  study,  a 
desired  berm  width  of  100  ft  from  the  construction  baseline  with  a  berm  height 
of  +5.4  ft  above  NGVD  was  applied  to  the  seven  profiles  representing  unique 
subreaches.  The  profile  closure  depth  of  -19  ft  NGVD  (Equation  5)  was  used  in 
these  calculations. 

33.  Using  ISRP,  a  volumetric  (ydVUnear  ft)  difference  was  computed 
between  the  representative  profiles  and  the  translated  profiles  (Table  7) 
(Appendix  D) .  This  volume  was  then  multiplied  by  the  length  of  subreach  rep¬ 
resented  by  that  particular  profile.  The  analysis  indicates  that  to  build  the 
23,090  ft  of  beach  to  design  dimensions  with  adjustment  to  depth  of  closure, 
1,228,000  yd^  of  borrow  material  with  a  grain  size  similar  to  the  native  beach 
would  be  required.  Theoretically,  when  redistribution  of  the  beach  fill  has 
occurred  across  the  active  profile,  forming  a  new  equilibrium  profile  shape. 

Table  7 

Design  Quantity  Estimates  for  Native  Beach  Subreach 


Subreach 

Length,  ft 

vd^/ft 

Yd^ 

Ra* .  vd^ 

5  th 

2,500 

8 

20,000 

31,000 

11th 

5,280 

56 

295,700 

458,300 

30th 

5,300 

57 

302,100 

468,200 

38th 

3,625 

81 

293,600 

455,100 

49  th 

1,975 

57 

112,600 

174,500 

55th 

1,850 

55 

101,800 

157,700 

61st 

40 

102.400 

158.700 

Total 

23,090 

1,228,200 

1,903,600 

*  Overfill  ratio  of  1.55  multiplied  by  cubic  yards  for  .each  subreach. 


the  entire  project  would  have  a  resulting  berm  width  of  approximately  100  ft 
from  the  baseline  with  a  berm  elevation  of  +5.4  ft  NGVD. 
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Offshore  borrow  material 

34.  Dredged  material  from  the  Thimble  Shoals  navigation  channel  was 
selected  by  CENAO  as  the  borrow  for  the  Virginia  Beach  Hurricane  Project  (Fig¬ 
ure  8) .  To  characterize  the  borrow  area,  24  cores  were  taken  within  Thimble 
Shoals.  Channel  sediment  samples  were  extracted  from  each  core  based  upon 
stratigraphic  units  or  unique  sediment  characteristics.  The  raw  grain  size 
data  from  each  unit  for  all  cores  was  provided  to  CERC  for  grain  size  com¬ 
posite  computations.  A  total  of  seven  zones  within  Thimble  Shoals  were  ini¬ 
tially  identified  by  CENAO  as  potentially  suitable  borrow  areas.  Table  8 

Table  8 

Cores  Samples  Used  to  Compute  Borrow  Area  Grain  Size  Composites 


Core  Number  Elevation,  ft  (NGVD) 


83VC-52 

-53.0 

to 

-60.0 

84VC-162 

-51.0 

to 

-53.3 

83VC-83 

-41.0 

to 

-45.0 

83VC-53 

-52.0 

to 

-55.4 

83VC-54 

-52.0 

to 

-55.5 

84VC-164 

-40.0 

to 

-54.0 

84VC-165 

-52.5 

to 

-54.5 

83VC-84 

-40.6 

to 

-55.5 

83VC-55 

-50.0 

to 

-54.0 

85VC-228 

-40.0 

to 

-55.0 

85VC-229 

-53.0 

to 

-59.0 

83VC-56 

-52.0 

to 

-60.0 

83VC-85 

-41.0 

to 

-55.0 

84VC-167 

-46.0 

to 

-56.0 

84VC-166 

-52.0 

to 

-56.0 

83VC-57 

-52.0 

to 

-60.0 

83VC-58 

-53.0 

to 

-60.0 

84VC-168 

-50.0 

to 

-58.0 

84VC-169 

-54.0 

to 

-61.0 

83VC-59 

-53.0 

to 

-60.0 

83VC-60 

-52.0 

to 

-60.0 

84VC-170 

-45.0 

to 

-53.0 

84VC-171 

-52.0 

to 

-57.0 

83VC-61 

-52.0 

to 

-57.0 

Borrow  Area 


#1 


#3 


#2  Y  j 

#4 

#5 


#6 


#7 


26 


r 


Figure  8.  Location  map  of  Thimble  Shoals  borrow  area 


indicates  cores  and  depth  locations  where  sediment  samples  were  taken  to  sta¬ 
tistically  create  composites  for  the  seven  zones. 
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35.  A  computer  program  was  developed  at  CERC  to  compute  composite  sedi¬ 
ment  grain  size  distributions  for  any  borrow  area  or  native  beach.  Various 
schemes  were  developed  to  weight  individual  sediment  samples  to  create  statis¬ 
tically  correct  composites.  Borrow  area  sediment  composites  could  be  computed 
in  this  program  based  upon  core  location,  depth,  beach  profile  location(s) ,  or 
any  combination  of  the  three.  For  this  particular  project,  borrow  area  com¬ 
posites  were  created  by  selecting  individual  samples  within  cores  between 
defined  depth  intervals.  The  percentage  of  s^lmple  within  the  specified  depth 
interval  was  multiplied  by  the  sample's  grain  size  distribution.  An  example 
would  be  to  create  a  borrow  area  composite  consisting  of  all  sediment  samples 
between  depths  of  -50  to  -60  ft  NGVD  taken  from  cores  83VC-52,  83VC-53,  and 
83VC-54.  If  a  particular  sample  in  one  of  the  cores  represented  a  strati¬ 
graphic  unit  from  -50  to  -55  ft,  a  weighting  factor  of  50  percent  was  multi¬ 
plied  by  the  Individual  weights  of  each  sieve  interval,  i.e.  (-50  -  -55)/ 

(-50  -  -60)  X  100  -  50  percent.  Remaining  samples  located  between  the  -50  to 
-60  ft  depth  interval  would  be  weighted  in  a  similar  manner.  Using  this  pro¬ 
cedure,  a  single  composite  was  created  for  each  core  with  each  core  composite 
possessing  a  weighting  factor  of  one.  The  result  was  a  single  grain  size  com¬ 
posite  which  represents  that  particular  borrow  area. 

36.  Table  9  represents  the  composite  grain  size  distribution  for  the 
four  borrow  areas  selected  by  CENAO.  Based  upon  results  from  the  seven  borrow 
sources,  four  borrow  sites  were  selected  for  further  investigation.  These 
sites  were  borrow  areas  designated  #2,  #3,  #6,  and  #7. 

Native  beach  material 

37.  Native  beach  sediment  characteristics  are  required  for  beach-fill 
design.  The  sediment  grain  size  characteristics  of  a  beach  vary:  (a)  across 
the  beach  profile,  (b)  along  the  beach,  and  (c)  between  seasons  (Hobson  1977). 
To  accurately  describe  this  variability,  temporal  and  spatial  native  beach 
sediment  samples  should  be  mathematically  combined  to  create  an  overall  com¬ 
posite  grain  size  distribution.  Using  this  procedure,  one  composite  grain 
size  distribution  reflects  the  three  components  of  variability.  Accuracy  of 
the  composite  is  crucial  as  the  native  beach  composite,  as  well  as  borrow 
composite,  determines  the  project's  overfill  ratio  and  ultimate  construction 
and  maintenance  costs. 

38.  Natural  native  sediment  characteristics  at  Virginia  Be<:  -h  have  been 
obscured  by  continuous  beach  replenishment  efforts  since  the  early  1950' s. 


28 


Subreaches 


Mean  SD* 


#2 

Mean  SD 
1.74  1.24 

Ra  R1 


1.72  1.27 

Ra  R1 


>ources _ 

#6 

Mean  SD 
1.78  1.26 

Ra  Ri_ 


1.11 

-Ei_ 


5  th 

March 

1.67 

0.86 

1.21 

0.63 

1.22 

0.59 

1.24 

0.64 

1.15 

0.73 

Sept 

1.56 

1.11 

1.16 

1.04 

1.15 

0.99 

1.19 

1.06 

1.15 

1.12 

Both 

1.64 

0.99 

1.16 

0.83 

1.16 

0.79 

1.19 

1.06 

1.15 

1.12 

11th 

March 

1.68 

0.84 

1.22 

0.60 

1.22 

0.55 

1.23 

0.60 

1.15 

0.69 

Sept 

1.49 

1.09 

1.23 

1.09 

1.21 

1.03 

1.26 

1.10 

1.23 

1.18 

Both 

1.62 

0.97 

1.18 

0.82 

1.18 

0.78 

1.21 

0.84 

1.11 

0.92 

30th 

March 

1.38 

0.90 

1.40 

0.95 

1.38 

0.89 

1.44 

0.97 

1.35 

1.08 

Sept 

1.63 

1.11 

1.11 

0.98 

1.11 

0,93 

1.14 

0.99 

1.07 

1.06 

Both 

1.50 

1.03 

1.23 

1.01 

1.22 

0.95 

1.27 

1.02 

1.20 

1.11 

38th 

March 

1.52 

0.85 

1.31 

0.74 

1.31 

0.68 

1.35 

0.75 

1.24 

0.86 

Sept 

1.44 

1.20 

1.32 

1.24 

1.27 

1.19 

1.36 

1.26 

1.37 

1.32 

Both 

1.46 

1.07 

1.26 

1.09 

1.24 

1.04 

1.30 

1.11 

1.27 

1.19 

49  th 

March 

1.47 

0.84 

1.36 

0.76 

1.35 

0.71 

1.39 

0.77 

1.29 

0.89 

Sept 

1.30 

1.09 

1.48 

1.29 

1.42 

1.23 

1.52 

1.31 

1.55 

1.40 

Both 

1.38 

0.98 

1.36 

1.07 

1.34 

1.01 

1.40 

1.09 

1.35 

1.19 

55th 

March 

1.46 

0.82 

1.38 

0.74 

1.37 

0.68 

1.41 

0.75 

1.31 

0.87 

Sept 

1.40 

1.04 

1.33 

1.12 

1.30 

1.06 

1.37 

1.14 

1.34 

1.23 

Both 

1.41 

0.94 

1.35 

0.98 

1.34 

0.92 

1.39 

1.00 

1.30 

1.10 

61st 

March 

1.45 

0.84 

1.37 

0.78 

1.36 

0.73 

1.41 

0.79 

1.31 

0.91 

Sept 

1.20 

1.19 

1.72 

1.51 

1.64 

1.44 

1.77 

1.53 

1.88 

1.61 

Both 

1.29 

1.06 

1.48 

1.27 

1.43 

1.21 

1.53 

1.29 

1.55 

1.39 

Winter 

1.59 

0.85 

1.27 

0.68 

1.26 

0.63 

1.30 

0.69 

1.20 

0.79 

Svimmer 

1.42 

1.13 

1.32 

1.20 

1.28 

1.14 

1.36 

1.22 

1.36 

1.29 

All 

1.55 

1.00 

1.30 

0.95 

1.27 

0.87 

1.34 

0.94 

1.25 

1.02 

SD  -  standard  deviation. 


Traverse  Group  Inc.  (1980)  performed  an  historical  sedimentological  study  at 
Virginia  Beach  using  data  from  1951  to  1977  and  found  several  trends.  The 

study  found  that  the  median  grain  size  above  mean  low  water  (MLW)  was  coarsest 

(0.38  mm,  1.40)  in  1951  gradually  becoming  finer  with  the  finest  (0.22  mm, 

2.19)  recorded  in  February  1977.  However,  the  median  grain  size  below  MLW  has 

remained  fairly  constant  through  time.  The  spatial  variation  for  the  five 
profiles  surveyed  is  Inconsistent  through  time.  The  study  also  found  the 
beach  sediment  has  become  more  uniform,  i.e.  lower  standard  deviation,  since 
1951 .  Sediment  characteristics  of  the  present  beach  reflect  mainly  the 
replenished  material.  These  trends  are  probably  a  result  of  the  continual 
beach  replenishment  practices. 

39.  To  describe  the  present  sedimentological  condition  at  Virginia 
Beach,  native  beach  sediment  data  for  13  profile  lines  were  provided  by  CENAO. 
Surface  sediment  samples  for  the  13  lines  were  collected  for  four  time  peri¬ 
ods:  March  1986,  September  1986,  March  1987,  and  September  1987.  Samples 
were  collected  at  eight  locations  along  each  profile  line:  +3.8,  1.8,  0.1, 
-1.6,  -6.6,  -11.6,  -16.6,  and  -21.6  ft  NGVD. 

40.  Native  composite  grain  size  distributions  were  computed  for  the 
seven  subreaches  identified  for  the  project.  These  seven  subreaches  were  the 
same  used  for  the  volumetric  analysis.  Unlike  the  profile  data  where  one 
profile  represents  a  subreach,  the  sediment  subreaches  were  represented  by 
data  from  three  profiles.  Sediment  data  from  three  profile  lines  were  used  to 
decrease  the  natural  variability  in  the  grain  size  statistics.  The  center 
profile  was  the  same  used  in  the  volumetric  analysis  with  additional  sediment 
data  provided  from  adjacent  profiles.  Winter  (March  1986  and  1987),  summer 
(September  1986  and  1987) ,  and  winter-summer  composites  were  computed.  The 
winter  and  summer  composites  represent  temporal  variations  of  the  native 
beach.  It  should  be  noted  that  winter- summer  composites  are  not  an  average  of 
the  two,  but,  a  composite  of  all  individual  samples  for  these  periods.  For 
native  beach  composites,  each  surface  sample  was  assigned  an  equal  weighting. 

41.  Typically,  sediment  samples  from  across  the  entire  profile,  i.e. 

+12  to  -30  NGVD,  are  included  in  the  grain  size  composite  for  the  native  mate¬ 
rial.  Stauble,  Hansen,  and  Blake  (1984)  found  that  the  inclusion  of  offshore 
samples  in  the  native  composite  tends  to  make  the  native  composite  finer, 
thereby  reducing  the  project's  overfill  ratio.  Subsequently,  less  volume  is 
apparently  needed  to  design  the  beach-to-project  dimensions.  Stauble,  Hansen, 
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and  Blake  suggest  using  only  the  intertidal  sediment  samples,  l.e.  between 
mean  high  and  low  tide,  to  create  the  native  beach  composite. 

42.  The  approach  of  Stauble,  Hansen,  and  Blake  (1984)  was  taken  on  this 
project  since  the  purpose  was  to  create  a  berm  that  could  withstand  extreme 
waves  and  run-up  generated  by  hurricanes  and  extratroplcal  storms.  By  using  a 
native  composite  composed  of  material  located  In  the  Intertidal  zone.  It  Is 
possible  to  characterize  the  particle  sizes  that  compose  the  active  wave 
rxm-up  zone  during  extreme  events.  The  chances  of  maintaining  integrity  of 
the  berm  are  Increased  If  these  particle  sizes  are  abundant  during  such  an 
event.  Therefore,  only  the  sediment  samples  from  +5.4,  3.4,  1.7,  0.0,  -5.0  ft 
MLU  were  used  to  create  the  native  beach  composite. 

43.  Table  9  represents  the  composite  mean,  standard  deviation  for  the 
seven  subreaches  on  the  native  beach.  Unlike  the  Traverse  Group  study,  a 
definite  spatial  gradation  In  the  mean  grain  size  at  Virginia  Beach  was 
observed  during  the  2  years  In  which  sediment  data  were  collected.  A  gradual 
Increase  In  the  composite  mean  diameter  was  Identified  from  south  to  north 
with  the  standard  deviation  remaining  fairly  constant.  Temporally,  the  March 
composites  tended  to  be  finer  that  the  September  composites. 

Beach-Fill  Design 

44.  Two  equations  are  typically  used  to  compute  losses  of  borrow  mate¬ 
rial  for  beach  nourishment  projects  (SPM  1984).  The  fill  factor  (Ra)  or  over¬ 
fill  ratio  estimates  the  volume  of  borrow  material  necessary  to  create  one 
unit  of  stable  beach  material  (James  1975).  For  example,  an  overfill  ratio  of 
2.0  Indicates  that  two  units  of  borrow  material  are  required  to  create  one 
unit  of  native  beach  material  after  winnowing  has  occurred.  A  second  equa¬ 
tion,  the  renourishment  factor  (Rj),  estimates  the  long-term  stability  of  a 
particular  borrow  relative  to  the  native  beach  material.  A  renourishment 
factor  of  two  suggests  that  the  borrow  material  Is  one -half  as  stable  as 
native  beach  material  or  that  renourishment  would  have  to  occur  twice  as  often 
to  maintain  specific  beach  dimensions.  This  factor  is  often  applied  to  select 
borrow  sources  with  similar  overfill  ratios  (Hobson  1977) , 

45.  As  previously  stated,  there  was  a  gradual  Increase  In  the  composite 
mean  grain  size  diameter  from  Rudee  Inlet  to  61st  Street.  Consequently,  over¬ 
fill  and  renourishment  factors  computed  for  each  subreach  reflect  this 
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gradation  by  increasing  from  south  to  north.  Table  9  represents  the  overfill 
and  renourishment  factors  for  all  subreaches  in  this  study. 

46.  The  overfill  factors  computed  for  each  subreach  could  have  been 
applied  separately  to  compute  project  volume  requirements.  However,  the  area 
of  most  concern  for  this  project  was  centered  in  the  vicinity  of  38th  to 
49th  Streets  where  there  has  been  historically  accelerated  erosion.  There¬ 
fore,  instead  of  applying  overfill  factors  for  each  subreach,  a  more  cautious 
approach  was  taken  by  applying  the  overfill  ratio  computed  for  the  area  of 
high  erosion  to  the  entire  project.  In  cooperation  with  CENAO,  an  overfill 
ratio  of  1.35  was  selected  for  the  entire  project  to  compensate  for  an 
expected  high  loss  of  finer,  unstable  material  during  the  initial  profile 
adjustment  period. 

47.  Computed  renourishment  factors  for  the  seven  subreaches  ranged  from 
0.63  to  1.69  (Table  9).  These  values  should  be  used  with  caution  since  the 
renourishment  equation  has  not  been  fully  verified.  Assessment  of  the  renour¬ 
ishment  factors  suggests  borrow  sites  #2,  #3,  and  #6  may  be  more  suitable 
borrow  sources  as  their  renourishment  factors  were  slightly  lower  compared 
with  borrow  site  #7.  An  economic  evaluation  based  upon  overfill  ratios, 
depth,  distance  to  beach,  and  available  volume  should  be  performed  prior  to 
selecting  the  final  borrow  source(s)  for  this  project. 

Design  volume  requirements 

48.  Due  to  the  physical  location  and  distance  of  Thimble  Shoals  borrow 
area  from  Virginia  Beach,  it  is  most  likely  that  borrow  material  will  be  exca¬ 
vated  using  a  hopper  dredge  and  then  pumped  out  onto  the  beach  via  an  offshore 
pumpout  station.  A  certain  percentage  of  borrow  material  is  lost  in  this  type 
operation  because  of  multiple  handling  (transfers)  of  the  material.  To 
account  for  this  anticipated  loss,  a  handling  loss  factor  (typically  10  to 

25  percent)  is  added  to  the  overfill  ratio  (Hobson  and  James  1979) .  In  a 
cooperative  decision  with  CENAO,  a  handling  loss  of  15  percent  was  applied  to 
the  project  overfill  ratio  of  1.35.  This  brings  the  final  overage  factor  for 
the  hurricane  protection  project  to  1.55.  Multiplying  this  factor  times  the 
volume  required  to  build  the  beach  to  design  dimensions  of  1,228,000  yd® 
yields  1,904,000  yd®  as  the  final  volume  of  borrow  material  required  to  obtain 
the  specified  design.  The  design  construction  profile  extends  205  ft  from  the 
construction  control  line  at  an  elevation  of  6.4  ft  NGVD  (8.0  MLW)  and  then 
tapers  to  existing  bottom  with  a  1:20  slope  (Appendix  D) . 
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Advance  nourishment  estimates 

49.  At  the  request  of  CENAO,  advanced  nourishment  volume  qtiantlties 
were  estimated.  Fortunately,  many  data  have  been  collected  at  Virginia  Beach 
concerning  historic  erosion  rates  and  subsequent  annual  nourishment  quanti¬ 
ties.  An  estimated  220,000  yd^  of  material  is  annually  bypassed  from  Rudee 
Inlet  and  placed  south  of  18th  Street.  North  of  18th  Street,  approximately 
150,000  yd^  of  material  is  annually  truck  hauled  and  deposited  on  the  beach. 
This  combined  quantity  of  370,000  ydVy®ar  tends  to  stabilize  a  majority  of 
beaches  for  a  period  of  approximately  1  year.  The  exception  is  an  area  in  the 
vicinity  of  38th  Street  that  tends  to  historically  erode  at  a  greater  rate 
than  the  adjacent  areas. 

50.  For  these  calculations,  it  was  assumed  the  design  volume  of 
1,904,000  yd®  would  provide  an  adjusted,  uniform  beach  width  of  100-ft  berm 
from  the  baseline  after  establishment  of  a  new  equilibrium  profile.  A  2 -year 
advance  nourishment  will  require  an  additional  740,000  yd®  (370,000  yd®  x 

2  years)  times  an  overfill  ratio  of  1.55  for  a  final  2 -year  advance  nourish¬ 
ment  quantity  of  1,145,000  yd®  (Table  10).  A  total  of  3,049,000  yd® 

(1,904,000  +  1,145,000  yd®^  of  material  would  be  required  for  the  entire  proj¬ 
ect  if  the  2 -year  advance  nourishment  is  performed  at  the  time  of  initial 
construction.  The  construction  profile  cross  sections  for  the  design  and  2- 
and  5 -year  advance  nourishment  quantities  are  represented  in  Appendix  E.  The 
2 -year  advance  nourishment  construction  profile  extends  290  ft  from  the  con¬ 
struction  control  line  and  then  tapers  to  existing  bottom  with  a  1:20  slope. 

An  overfill  ratio  was  applied  to  the  advance  nourishment  quantity  since  it  was 
assumed  that  losses  similar  to  the  design  volume  would  occur.  The  2-year 
advance  nourishment  will  provide  approximately  an  additional  32  yd®/linear  ft 
of  beach  after  winnowing  of  fines  and  handling  losses.  This  volume  of  mate¬ 
rial  can  be  roughly  translated  into  32  ft  of  addition  berm  width  using  the 
rule  of  thumb  where  1  yd®  equals  1  ft  of  beach  width. 

51.  Similar  calculations  were  performed  for  a  5-year  advance  nourish¬ 
ment.  A  5 -year  advance  nourishment  would  require  approximately  1,850,000  yd® 
(370,000  yd®  x  5  years)  times  an  overfill  ratio  of  1.55  for  a  5-year  advance 
nourishment  quantity  of  2,867,000  yd®  (Table  11)  (Appendix  E) .  A  total  of 
4,777,000  yd®  of  (1,904,000  +  2,867,000  yd®)  material  would  be  required  for 
the  entire  project  if  the  5 -year  advance  nourishment  is  performed  at  the  time 
of  initial  construction.  This  volume  of  material  can  be  translated  into  an 
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Table  10 


Two-Year  Advance  Nourishment  Quantity  Estimates 
for  Native  Beach  Subreach 


Subreach 

Length 

ft 

vd^/ft 

vd^ 

Ra.*  vd^ 

Design  & 
Advance 

vd^ 

5  th 

2,500 

32 

80,000 

124,000 

155,000 

11th 

5,280 

32 

169,000 

261,900 

720,200 

30th 

5,300 

32 

169,600 

262,900 

731,100 

38th 

3,625 

32 

116,000 

179,800 

634,900 

49  th 

1,975 

32 

63,200 

98,000 

272,500 

55th 

1,850 

32 

59,200 

91,800 

249,500 

61th 

2.560 

32 

82.000 

127.000 

285.700 

Total 

23,090 

739,000 

1,145,000 

3,048,900 

*  Overfill  ratio  of  1.55  multiplied  by  cubic  yards  for  each  subreach. 


Table  11 

Five-Year  Advance  Nourishment  Quantity  Estimates 
for  Native  Beach  Subreach 


Design  & 


Subreach 

Length 

ft 

vd^/ft 

vd^ 

Ra.*  vd^ 

Advance 

vd^ 

5  th 

2,500 

80.1 

200,300 

310,400 

341,400 

11th 

5,280 

80.1 

423,000 

655,500 

1,113,800 

30th 

5,300 

80.1 

424,600 

658,000 

1,126,200 

38th 

3,625 

80.1 

290,400 

450,000 

905,100 

49th 

1,975 

80.1 

158,200 

245 , 200 

419,700 

55th 

1,850 

80.1 

148,200 

229,700 

387,400 

61th 

2. 5.^9 

80.1 

205.100 

317.800 

476.500 

Total 

23,090 

1,849,800 

2,866,600 

4,770,100 

*  Overfill 

ratio  of  1 

.55  multiplied  by  cubic  yards 

for  each  subreach. 
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additional  80  ft  of  berm  width  provided  by  this  quantity.  The  5 -year  advance 
nourishment  construction  profile  extends  410  ft  from  the  construction  control 
line  and  then  tapers  to  existing  bottom  with  a  1:20  slope. 


PART  IV;  SUMMARY  AND  CONCLUSIONS 


52.  Model  simulations  show  that  the  proposed  100-ft  berm  of  the  design 
dune /beach  profile  for  the  no -seawall  case  is  adequate  for  protection  from 
either  a  single  or  multiple  design  storm  event.  Additional  model  results 
indicate  that  a  minimum  width  of  50.0  ft  would  provide  protection  against 
breaching  by  the  design  storm  although  erosion  of  the  crest  of  the  dune  would 
occur.  This  analysis  assumes  that  the  design  berm  configuration  is  maintained 
at  all  times  and  that  the  1-D  assumption  is  not  severely  violated.  If  the 
verification  results  reported  by  Birkemeier  et  al.  (1987)  are  taken  into  con¬ 
sideration,  a  factor  of  safety  of  approximately  2.0  should  be  applied  to  the 
present  analysis.  This  figure  results  from  the  fact  that  erosion  in  5  of  the 
13  cases  was  underestimated  by  55  percent  or  that  the  observed  erosion  volume 
was  1.82  times  the  predicted.  A  factor  of  safety  of  2.0  was  also  recommended 
by  Kraus  et  al.  (in  preparation)  in  a  dune  erosion  model  application  to  the 
shoreline  of  New  Jersey.  In  view  of  these  results  and  the  fact  that  wide¬ 
spread  damage  may  occur  if  the  dune  is  breached,  the  100-ft  berm  width  pro¬ 
posed  by  CENAO  is  recommended.  According  to  the  stage -frequency  curve  computed 
by  CENAO,  this  would  provide  protection  for  the  100-year  storm  (i.e.  the  August 
1933  storm)  and  provide  a  reasonable  margin  of  safety. 

53.  The  proposed  dune/beach  profile  with  a  100-ft  berm  width  for  the 
design  case  with  a  seawall  also  provides  adequate  protection  against  the  design 
storm.  For  the  seawalled  case,  the  berm  will  cease  receding  when  recession 
reaches  the  seawall.  At  this  point,  scour  at  the  face  of  the  seawall  will 
begin.  Additional  simulations  were  performed  in  which  the  berm  width  was 
reduced  to  50.0  ft  in  order  to  evaluate  the  potential  erosion  in  front  of  the 
seawall.  Results  of  the  multiple  design  storm  simulation  showed  that  the  berm 
elevation  was  lowered  over  its  entire  width.  However,  the  elevation  loss  at 
the  seawall  was  not  sufficient  to  endanger  seawall  integrity  and  was  comparable 
to  the  value  of  approximately  2  ft  of  erosion  used  in  computing  overtopping 
quantities  in  the  physical  model  study  (Lillycrop  et  al.  1988).  In  view  of 
this  result  and  the  variability  factor,  the  100-ft  width  at  elevation  +5.4  NGVD 
design  is  recommended  as  an  effective  measure  of  protection. 

54.  The  profile  translation  method  was  used  to  determine  the  quantity 
of  material  needed  to  build  the  beach,  in  its  present  condition,  to  a  beach 
possessing  a  berm  width  of  100  ft  after  the  establishment  of  a  new  equilibrium 
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profile.  The  CERC  was  provided  with  four  sets  of  profile  data  spanning  a 
period  of  2  years.  The  March  1987  profile  data  set  was  selected  to  calculate 
the  design  volume  for  the  project.  This  data  set  provided  the  most  conserva¬ 
tive  (greatest)  volumetric  results.  Seven  unique  subreaches  of  the  project 
beach  were  identified  and  represented  by  seven  profile  lines.  Total  volume 
necessary  to  construct  the  design  berm  without  advance  nourishment  was  calcu¬ 
lated  to  be  1,227,000  yd^. 

55.  Composite  grain  size  analysis  was  conducted  for  seven  borrow  areas 
within  Thimble  Shoals,  and  four  areas  were  selected  for  further  Investigation. 
In  a  similar  manner,  composite  grain  size  analysis  were  computed  for  the  same 
seven  native  beach  subreaches  used  in  the  volumetric  assessment.  Overfill 
ratios  and  renourishment  factors  were  calculated  using  all  native  beach  and 
borrow  area  composites.  Overfill  ratios  ranged  from  1.07  to  1.88  for  the 
various  scenarios  of  native  and  borrow  material.  Since  the  area  of  acceler¬ 
ated  erosion  is  in  the  vicinity  of  38th  Street,  an  overfill  ratio  of  1.35  was 
selected  by  CERC  and  CENAO  to  be  applied  for  the  entire  project.  Because  of 
the  type  of  dredging  operation  anticipated  for  Thimble  Shoals,  a  handling  loss 
of  15  percent  was  applied  to  the  overfill  ratio  yielding  a  final  overage  fac¬ 
tor  of  1.55  for  the  entire  project.  Multiplying  the  final  overage  factor 
(1.55)  times  the  volume  of  material  required  to  build  the  beach  to  design 
dimensions  (1,227,000  yd^)  yields  1,902,000  yd®  of  material  to  build  the  proj¬ 
ect  to  specifications.  Evaluating  the  renourishment  factors  suggests  borrow 
areas  #2,  #3,  and  #6  are  probably  the  most  suitable;  however,  further  economic 
evaluation  should  be  performed  before  making  the  final  borrow  site  selection. 

56.  Currently,  an  estimated  370,000  yd®/year  of  material  is  artifi¬ 
cially  placed  in  the  project  area  on  a  yearly  basis  to  stabilize  the  beach 
width.  Based  upon  current  stabilization  practices,  2-  and  5-year  advance 
nourishment  quantities  were  estimated.  A  2-  and  5-year  advance  nourishment 
would  require  an  additional  1,147,000  and  2,867,500  yd®  of  material,  respec¬ 
tively,  to  the  design  volume  of  1,902,000  yd®  giving  a  total  of  3,049,000  and 
4,769,500  yd®,  respectively.  These  quantities  would  theoretically  stabilize 
the  beach  at  100  ft  from  the  construction  control  line  for  the  respective  time 
interval.  An  overfill  ratio  of  1.55  is  included  in  the  advance  nourishment. 
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APPENDIX  A:  STOKM  SURGE  CURVES 

(Storm  surge  curves  used  in  the  numerical  modeling  test  represent  three 
northeasters  and  two  hurricanes . ) 
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APPENDIX  B:  PRESTORM  AND  POSTSTORM  PROFILES 

(Prestorm  and  poststorm  profiles  are  presented  with  and  without  a  protective 
seawall.  Elevations  are  relative  to  National  Geodetic  Vertical  Datum.) 
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APPENDIX  C:  REPRESENTATIVE  PROFILES  FOR  SUBREACHES 

(Representative  profiles  for  the  seven  native  beach  subreaches  are  presented. 
Horizontal  distance  is  relative  to  the  construction  control  line.  Elevations 
are  relative  to  mean  low  water  (MLW) . ) 
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APPENDIX  D:  DESIGN  AND  DESIGN  CONSTRUCTION  PROFILES 

(Presented  are  the  design  profiles  and  design  construction  profiles  for  the 
seven  native  beach  subreaches.  Horizontal  distance  is  relative  to  the 
construction  control  line.  Elevations  are  relative  to  mean  low  water  (MLW) . ) 
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APPENDIX  E!  ADVANCE  NOURISHMENT  PROFILES 

(Design  construction  profiles  and  2-  and  5 -year  advance  nourishment  profiles 
are  presented  for  the  seven  native  beach  subreaches.  Horizontal  distance  is 
relative  to  the  construction  control  line.  Elevations  are  relative  to  mean 
low  water  (MLW) . ) 
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